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RESUMEN 
ABSTRACT 

The age distribution of the central stars of planetary nebulae (CSPN) is es- 
timated using two methods based on their kinematic properties. First, the 
expected rotation velocities of the nebulae at their Galactocentric distances 
are compared with the predicted values for the rotation curve, and the dif- 
ferences are attributed to the different ages of the evolved stars. Adopting 
the relation between the ages and the velocity dispersions determined by the 
Geneva-Copenhagen survey, the age distribution can be derived. Second, the 
U, V, W, velocity components of the stars are determined, and the corre- 
sponding age- velocity dispersion relations are used to infer the age distribu- 
tion. These methods have been applied to two samples of PN in the Galaxy. 
The results are similar for both samples, and show that the age distribution of 
the PN central stars concentrates in ages lower than 5 Gyr, peaking at about 
1 to 3 Gyr. 
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1. INTRODUCTION 

Planetary nebulae (PN) are evolved objects ejected by stars with main 
sequence masses in the range of 0.8 and 8 M©, so that the expected ages of 
their central stars are of the order or greater than about 1 Gyr. However, the 
relatively large mass bracket of their progenitor stars implies that an age dis- 
tribution is to be expected, which has some consequences in the interpretation 
of the PN data in the Galaxy and other stellar systems. The determination 
of ages of the central stars is a difficult problem, and most usual methods 
have large uncertainties when applied to intermediate and old age objects. 
We have recently developed three different meth ods t o estimate the age dis- 
tribu t ion of the CSPN (Maciel, Costa & Idiart boiol see also Maciel et al. 



20031 l2005l I2006I ). and applied these methods to a sample of PN in the disk 
of the Galaxy, most of which are located in the solar neighbourhood, within 
3 kpc from the Sun. These methods include the determination of the age 
distribution of CSPN using (i) an age-metallicity relation that also depends 
on the Galactocentric distance, (ii) an age-metallicity relation obtained for 
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the disk, and (iii) the central star masses obtained from the observed nitrogen 
abundances. We concluded that most CSPN in our sample have ages under 
6 Gyr, and that the age distribution is peaked around 2-4 Gyr. The aver- 
age uncertainties were estimated as 1-2 Gyr, and the results were compared 
with the expected distribution based both on the observed mass distribution 
of white dwarfs and on the age distribution derived from available masses of 
CSPN. 

In the present work we develop two additional and more accurate methods 
to estimate the age distribution of the CSPN based on their kinematical prop- 
erties, namely: (i) A method based on the expected rotation velocities of the 
nebulae at their Galactocentric distances, which are then compared with the 
predicted values for a given rotation curve; the differences are attributed to 
the different ages of the evolved stars; (ii) A method based on the derived U, 
V, W, velocity components of the stars and their corresponding dispersions. In 
both cases, the age- velocity dispersion relations from the Geneva-Copenhagen 
survey are used to infer the age distribution. These methods are applied to 
two P N sani ples, (i) the previous sample of disk PN used by Maciel, Costa & 
Idiart (l2010fl . for which a detailed data set is available, and (ii) a sample con- 



taining all PN for which accurate radial velocities are known. The methods 
are developed in Section 2, and the samples used are described in Section 3. 
The main results and discussion are given in Section 4. 

2. DETERMINATION OF THE AGE DISTRIBUTION OF CSPN 

2.1. Method 1: The PN rotation velocity 

As objects of intermediate age, PN the disk of the Galaxy describe a rota- 
tion curve similar to the one defined by younger objects, such as HII regions, 
althou gh with a higher dispersion, as discussed in detail by Maciel and Lago 



(|2005f ). Therefore, the discrepancies between the rotation velocities inferred 
from the PN radial velocities and distances and the expected velocities from 
the known rotation curve may be at least partially ascribed to their evolved 
status. In other words, a given nebulae located at a distance d, with galac- 
tic coordinates I and h and observed heliocentric radial velocity Vr{hel) can 
be associated with a rotation velocity 0{R), after obtaining its Galactocen- 
tric distance R and its radial velocity relative to the Local Standard of Rest 
(LSR), Vr{LSR). Assuming circular orbits, the rotation velocity 0{R) at the 
Galactocentric distance R can be written as 



OiR) = ^ 

Ho 



Vr{LSR) 

sini cos 6 



(1) 



where Rq and ^re the Galactocentric distance and rotation velocity at th e 
solar position (see for example Maciel & Lago ^2005, Maciel & Dutra 1992f) . 



On the other hand, the expected rotation velocity at the given Galactocentric 
distance, dc{R), can be obtained from an adopted rotation curve. The differ- 
ence A9 = l^'(^) — Sc{R)\ can then be considered as proportional to the age 
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COEFFICIENTS OF THE POLYNOMIAL GIVEN BY EQ. (3). 
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R (kpc) - 0.765 0.765 - 2.9 2.9 - 3.825 3.825 - 13 > 13 



ao 


0.0 


325.0912 


329.8 


-2346.0 


ai 


3069.81 


-248.1467 


-250.1 


2507.60391 


ai 


-15809.8 


231.87099 


231.87099 


-1024.068760 


as 


43980.1 


-110.73531 


-110.73531 


224.562732 


a4. 


-68287.3 


25.073006 


25.073006 


-28.4080026 


as 


54904.0 


-2.110625 


-2.110625 


2.0697271 


ae 


-17731.0 






-0.080508084 


aj 








0.00129348 



230.6 



difference between the PN and the objects defining the rotation cur ve. W e 



1998 ) 



20081) . 



have adopted the radial velocities from t he cata logue by Durand et al. 
and two distance scales, those by Maciel (|l984l ) and Stanghellini et al. 
The first one was based on a relationship between the ionized mass and the 
radius of th e nebu lae, while the second is an update of the distance s cale b y 
Cahn et al. (jl992l ) , using a modified Shklovksy method following Daub (jl982i ) . 
Since the distances of planetary nebulae in the Galaxy may contain large in- 
dividual uncert aintie s, the use of two different scal es whi ch are considered as 
"short" (Maciel 1984 ) and "long" (Stanghellini et al. 20081) warrants that these 
uncertainties will not affect the derived age distributions. We have adopted 
i?o = 8.0 kpc for the distance of the Sun to the centre and = 220km/s for 
the rotation velocity at Rq. Sligh tly diff erent valu es can be found in the liter- 
ature (see for example Ferryman, 20091 and Reid, 2O10l) . but the values above 
are frequently adopted, so that a comparison with other works is made easier. 
For the "theoretical" rotation curve we have also adop ted two possibilities, 
namely, the PN derived curve by Maciel & Lago |2005l ). and the HII region 
derived curve by Clemens |l985l) . In the first case, the rotation velocity can 
be written as 



0c{R) = ao + ai i? + 02 i?^ 



(2) 



where the constants are oq = 269.2549, ai = —14.7321, and 02 = 0.7847, the 
Galactocentric distance R is in kpc and 9c{R) in km/s. For the CO/HII region 
based Clemens (|l985l ) curve, we have made an adjustment for Rq = 8.0 kpc 
and Oo{R) = 220 km/s, in which case we have 



0,{R) =Y,a, W . (3) 
where the constants are given in Table 1, with the same units as in Eq.(2). 
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TABLE 2 
COEFFICIENTS OF EQ. (4). 





a 


b 


u 


0.39 


1.31 


V 


0.40 


1.10 


w 


0.53 


0.94 


Total 


0.40 


1.40 



The recent Geneva-Copenhagen S urvey of the Solar Neighbourhood (cf. 



Nordstrom et al. l2004i Holmberg et al. I2007ll2009l ) has considerably improved 
the relations involving the ages, kinematics, and chemical composition of a 
large sample containing about 14000 F and G nearby stars. Using basically 
the original Hipparcos parallaxes, uvby/3 photometry and the Padova stellar 
evolution models, several basic relations were investigated. In particular, high 
correlations have been obtained between the velocity dispersions au, ay, crw, 
and ax and the age of the star, which clearly show a smooth increase of the 
velocity dispersions in the U, V, W componen ts and total velocity T with 
time. From the calibration by Holmberg et al. (|2009l ) these correlations can 
be approximately written as 



log cr = a log t + 6 , (4) 

where the age t is in Gyr and the constants a, b are given in Table 2. This 
approximation is valid in the age interval < i(Gyr) < 14 with an estimated 
average age uncertainty of about 25%. Method 1 consists of assuming that 
the discrepancy in the rotation velocity A9 is due to the evolved status of the 
CSPN, so that we should expect a correlation between A9 and the velocity 
dispersion, as given by Eq. (4). Since in this method we are using the rotation 
velocity, we have considered two possibilities, according to which the velocity 
discrepancy A6 can be associated with (i) the V component of the total ve- 
locity (cry), or (ii) the total velocity (ctt)- Moreover, since we are adopting 
two distance scales and two theoretical rotation curves, we have 8 different 
age distributions for Method 1, characterized by the timescales ti to ts, as 
explained in Table 3. 



2.2. Method 2: The U, V, W velocity components 

Method 2 is also a kinematic method, and in principle more accurate than 
Method 1, as discussed in more detail in Section 4. From the PN radial 
velocities and distances, we have estimated their proper motions both in right 
ascension and declination, /i^ and fig. We have assumed that, in average, the 
tangential velocities are similar to the radial velocities, namely Vt — Vr- In 
view of the large distances of the nebulae, this hypothesis in practice does 
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TABLE 3 
PARAMETERS FOR METHOD 1. 



Distance 


Rotation Curve 


Dispersion 


Age 


Maciel 


PN 


av 


ti 


Maciel 


PN 


(7T 


ti 


Maciel 


Clemens 




ts 


Maciel 


Clemens 


(JT 


U 


Stanghellini 


PN 


av 


U 


Stanghellini 


PN 


(7t 


te 


Stanghellini 


Clemens 


ov 


h 


Stanghellini 


Clemens 







not introduce any major uncertainties in the results. Considering further 
the equatorial coordinat es (a, 6) of the PN, we have used the equations by 



Boesgaard and Tripicco (|l986l ) to derive the U, V, W velocity components 
of the nebulae, as well as the total velocity T and the velocity dispersions 
ajj, cry, cFw, and <jt- According to these equations we derive the following 
parameters: C = /(d), X = f{C,^ia,Hs,a,S,Vr), Y = /(C, ^q, /^a, a, 5, K), 
and Z — f{C,fis,S,Vr), from which the velocities can be written as U = 
f{X,Y,Z), V = f{X,Y,Z), W = f{X,Y,Z), and T = f{X,Y,Z), so that the 
dispersions are given by 



^/i^^'^' (5) 



where Vi stands for the velocities U, V, W, T. Then, we have again used the 
detailed correlations between the velocity dispersio ns and the ages as given by 



the Geneva-Copenhagen survey (Holmberg et al. 120091 ) . adopting the same 



coeffi cients given in Table 2. We ha ve used the same distance scales (Maciel 



19841 and Stanghellini et al. 120081 ) . so that we have again 8 different age 
distributions, corresponding to the timescales tg to ^ig, as described in Table 4. 

In practice, we have considered several additional cases, in order to better 
investigate the hypothesis of Vt ~ Vr- Assuming that these velocities are of 
the same magnitude, but allowing for the possibility of different signs, we 
have as a result several possibilities for the proper motions fXa and fis, all of 
which are consistent with either Vt ~ or \Vt \ ~ \ Vr\. It turns out that these 
possibilities produce very similar age distributions, which will be discussed in 
Section 4. Therefore, we will present only the distributions of the ages tg to 
ii6, as defined in Table 4, for the cases where c^i 0. 

An interesting alternative to overcome the lack of proper motion and tan- 
gential velocity measurements would be to ap ply th e singular value decompo- 
sition (SVD) technique, as used by Branham (|2010l ) to solve the inverse prob- 
lem, that is, obtaining the space velocities from available proper motions. 
However, in view of the similarity of the results for different assumptions 
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TABLE 4 
PARAMETERS FOR METHOD 2. 



Distance Dispersion Age 



Maciel 




tg 


Maciel 


ay 


tio 


Maciel 


aw 


tu 


Maciel 


Gt 


tl2 


Stanghellini 


ou 


tl3 


Stanghellini 


av 


tl4 


Stanghellini 


aw 


tl5 


Stanghellini 


aT 





regarding the tangential velocities, it is unlikely that this technique would 
produce very different results than presented here. 

3. THE SAMPLES 



As mentioned in the Introduction, we have considered two samples of Milky 
Way PN. In order to make comparisons with our pr evious work , we h ave first 
considered the same sample used by Maciel et al. ( 2003 , 20051 2006 ). which 
we will call Sample 1. This sample contains 234 well-observed nebulae located 
in the solar neighbourhood and in the disk, for which all data were obtained 
with the highest accuracy. Their Galactocentric distances are in the range 
4 < i?(kpc) < 14, and most (69%) are located in the solar neighbourhood, 
with distances d < 3kpc. 

The second sample considered in this work, called Sample 2, includes all 
the nebulae for whi ch ac curate radial velocities are available in the catalogue 
by Durand et al. (|l998l ). comprising 867 objects. This is a more complete 
sample, so that it is expected that the derived results can be extended to 
the observed population of PN in the Galaxy. In both samples, the number 
of nebulae used depends on the availabil ity of the statistical distances . The 
actual number of objects from the Maciel (1984) and Stanghellini et al. ( 20081 ) 
distance scales are 195 and 170 for Sample 1 and 493 and 403 for Sample 2, 
respectively. We have then applied the approximation given by Eq. (4) for 
both samples, with the coeSicients shown in Table 2, considering only the 
objects for which ages in the interval < i(Gyr) < 14 could be obtained. 



4. RESULTS AND DISCUSSION 



The main results for the age distribution of the CSPN are shown in Figures 
1-4, where we have used the age parameter definitions given in Tables 3 and 
4 for Methods 1 and 2, respectively. Figures 1 and 2 refer to Sample 1, while 
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TABLE 5 

FRACTION OF STARS AT THREE AGE INTERVALS. 





At (Gyt) 


0-3 


3-6 


> 6 


Method 1 


ti 


0.57 


0.13 


0.30 






0.62 


0.18 


0.20 




ta 


0.57 


0.19 


0.24 




t4 


0.67 


0.18 


0.16 




ts 


0.51 


0.13 


0.36 




U 


0.71 


0.17 


0.12 




tr 


0.61 


0.15 


0.24 




ta 


0.71 


0.11 


0.18 


Method 2 


tg 


0.76 


0.12 


0.12 




tio 


0.79 


0.10 


0.11 




til 


0.92 


0.04 


0.04 




tl2 


0.77 


0.18 


0.05 




tis 


0.78 


0.10 


0.12 




tl4 


0.78 


0.11 


0.11 




tl5 


0.93 


0.03 


0.04 




tie 


0.76 


0.18 


0.06 



figures 3 and 4 refer to Sample 2. It can be seen that the age distributions 
obtained by both methods are similar, in the sense that most objects have 
ages under 5 Gyr, with a strong peak at ages typically between 1 and 3 Gyr. 
The histograms of Figures 3-4 are summarized in Table 5, where the fraction 
of stars obtained by Method 1 (ages ti to ts) and Method 2 (ages tg to tie) 
are shown for three age bins, namely — 3 Gyr, 3 — 6 Gyr, and t > 6 Gyr. 

The similarity of the results of both methods is remarkable, especially con- 
sidering that Method 2 is probably more accurate than Method 1. Method 2 
consists of straightforward calculations of the velocities and velocity disper- 
sions followed by an application of relatively accurate correlations involv- 
ing the kinematics and ages of the objects considered. On the other hand. 
Method 1 is based on the assumption that the differences between the observed 
and predicted rotation velocities are essentially due to age effects. However, 
other processes may be important, such as deviations from the circular rota- 
tion, which is particularly important for nearby objects. According to Table 5, 
in all cases the vast majority of CSPN have ages under 3 Gyr. For Method 1 
the total fraction of objects with t < 3 Gyr is 50 — 70%, while for Method 2 
this fraction is somewhat higher, 70 — 90%. It is unlikely that this is a result 
from biased samples, as the results for the larger Sample 2 are essentially the 
same as in the smaller Sample 1. It should be pointed out that the latter, 
albeit smaller, includes only well studied nebulae, for which all individual 
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parameters (distances, velocities, abundances) are better determined. 

Also, there are no significant differences in the results using the different 
velocity components U, V, W, and T. For Method 1, the distributions using 
the V velocity component are essentially the same as those using the total 
velocity, for both distance scales and samples. For Method 2, the distributions 
are slightly more concentrated in the first few age bins for the W component, 
compared with the distributions for the U and V components and the total 
velocity, again for both distance scales and samples. Since the W component 
is more clearly associated with the disk heating, essentially caused by age 
effects, the corresponding distributions are probably more accurate. 

Similar remarks can be made regarding the adopted values for the proper 
motions. As mentioned at the end of Section 2, the results shown here assume 
that fia — IJ-s — 0. Adopting nonzero values for these quantities {fia — 
Us 7^ 0), either the V oi W component distributions become slightly less 
concentrated at the first few age bins, but most objects still have ages under 
about 4 Gyr. Again, the application of the SVD technique could be useful to 
confirm these results. 

The uncertainties in the distances of the Milky Way FN are difficult to 
estimate, but the procedure adopted here ensures that the obtained age distri- 
butions are not particularly affected by the individual distances of the objects 
in the samples. As mentioned in Section 2, we have adopted two very differ- 
ent statistical scales, and the derived age distrbutions are essentially the same 
in both cases. The individual distances may depend on the particular scale, 
but the results shown in Figures 1-4 and in Table 5 do not depend on the 
choice of the distance scale. This can be seen by comparing the results for the 
timescales ti — with those for — ts, or the results for ig — ti2 with those 
for ti3 - tie. 

The uncertainties in the radial velocities also do not seem to have an 
i mpor tant effect on the age distributions. In the catalogue by Durand et al. 
(|l998l ). most objects (~ 90%) have uncertainties smaller than 20 km/s, and 
many objects h ave n iuch lower uncertainties. Concerning Method 1, from 
Maciel & Lago (|2005l) . the average rms deviation in the rotation velocity is 
about 50 km/s for FN, which can be c omp ared with the values of abo ut 20 
km/s for HII regions (see also Clemens Il985l and Maciel & Dutra 119921 ). 

Frobably the main uncertainty of the age distributions is due to the cali- 
bration between the stellar ages and the velocity dispersions, given by Eq. (4), 
which affects both Method 1 and 2. From the Geneva- Copenhagen Survey, 
this relation has a dispersion of about 20 km/s in average, which corresponds 
roughly to an age uncertainty of about 25%, amounting to less than 1.2 Gyr 
for the objects of Figures 1-4. Therefore, the uncertainties of the present 
method are comparable and probably smaller than in the case of the methods 
based on age-metallicity relations considered by Maciel et al. (l2010f ). 

The results for Sample 2 are not essentially different from those of Sam- 
ple 1, so that a direct comparison can be made with the results by Maciel et al. 
(|2010l ). The results of both investigations are similar, even though the present 
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methods arc com pletel y independent of the metallicity-based methods used 
by Maciel et al. |2010l ). The main difference is that the kinematic methods 
used in the present investigation suggest somewhat lower ages for the CSPN 
in our samples. In this respect, these results fit nicely to the probability dis- 
tribution for the progenitors of the CSPN according to Maciel et al. ( 20ld cf. 
figure 7, dashed line). In this case the well known relatio n bet ween the main 
sequence mass and the stellar ages by Bahcall & Piran (jl983l ) was adopted, 
taking t — 10 Gyr for IMq stars on the main sequence. Taking into account 
the uncertainties of the methods, which are typically in the range 1 — 2 Gyr, 
this case was considered as the most realistic, so that it is reassuring that the 
kinematic methods produce similar results. 
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Fig. 1. Age distribution of CSPN, Method 1, Sample 1. 



KINEMATIC AGES OF CSPN 



11 






age (Gyr) age (Gyr) 

Fig. 2. Age distribution of CSPN, Method 2, Sample 1. 
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Fig. 3. Age distribution of CSPN, Method 1, Sample 2. 
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Fig. 4. Age distribution of CSPN, Method 2, Sample 2. 



